Purifying GHZ States Using Degenerate Quantum Codes 
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Degenerate quantum codes are codes that do not reveal the complete error syndrome. 
Their ability to conceal the complete error syndrome makes them powerful resources in 



W). 

certain quantum information processing tasks. In particular, the most error-tolerant way to 



purify depolarized Bell states using one-way communication known to date involves degen- 
erate quantum codes. Here we study three closely related purification schemes for depolar- 



ized GHZ states shared among m > 3 players by means of degenerate quantum codes and 



one-way classical communications. We find that our schemes tolerate more noise than all 
other one-way schemes known to date, further demonstrating the effectiveness of degenerate 
quantum codes in quantum information processing. 
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Quantum error correcting codes, unlike their classical counterparts, may not reveal the complete 

HQ. 



error syndrome. Codes with this property are known as degenerate codes HI, |2|]. In a sense, 
degenerate codes pack more information than non-degenerate ones because different quantum 
errors may not take the code space to orthogonal spaces. By carefully utilizing the degenerate 
property, degenerate codes are useful resources in quantum information processing. Examples 
showing their usefulness were provided by Shor and his co-workers Q]. In particular, they 
showed that a carefully constructed degenerate code is able to purify Bell states passing through 
a depolarizing channel with fidelity greater than 0.80944 [2]. Their scheme is more error-tolerant 
than all the known one-way depolarized Bell state purification schemes involving non-degenerate 
codes to date. 

It is instructive to ask if the degenerate codes can be used to improve the error-tolerant level 
of existing one-way multipartite purification protocols. Here we provide such an example by 
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considering the purification of shared GHZ states. Specifically, suppose that a player prepares 
many copies of perfect GHZ state in the form 

IO" ,+ ) = — (|0® m > + 11®")) . (1) 
V2 V ; 

For each perfect GHZ state, he/she keeps one of the qubit and sends the other to the remaining 
players through a depolarizing channel so that upon reception of their qubits, these m players share 
copies of Werner state 



W F = F |O m+ ) <O m+ | + (/ - |O m+ ) <O m+ |) , (2) 

where F is the fidelity of the channel and / is the identity operator. Now, the players wanted to 
distill shared perfect GHZ states using an one-way purification scheme that works for as small 
a channel fidelity as possible. Clearly, this task is a generalization of the Bell state distillation 
problem investigated by Shor and his co-workers yj, [2Q. 

We begin our study by defining a few notations and reviewing prior arts in Sec. [Til Then we 
introduce three closely related one-way multipartite purification protocols involving concatenated 
degenerate codes and analyze their performances in Sec. Hill Actually, all three protocols use the 
same repetition code as their inner codes. Moreover, in the case of m = 2, one of the our protocols 
is a generalization of the scheme proposed by DiVincenzo et al. [2]. Most importantly, for m > 3, 
our protocols are the most error tolerant ones discovered so far in the sense that ours can distill 
shared GHZ states from copies of Werner state in the form of Eq. © with a fidelity F so low that 
no other one-way purification schemes known to date can. Our schemes can also be generalized to 
the case when the Hilbert space dimension of each quantum particle is greater than 2. We briefly 
discuss this issue in Sec. [TV] Finally, we summarize our findings in Sec.lVl 

II. PRIOR ARTS 
A. Some notations 

Given that m > 2 players share N noisy GHZ states in the form of Eq. (0Q). Clearly, the GHZ 
state is stabilized by its stabilizer generators, namely, 

S o = XqX\ ■ ■ ■ X m -\ , 

Si = Z Zi (3) 
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for 1 < i < m - 1, where 
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denote the spin flip and phase shift operation acting on the ith qubit respectively. For simplicity, 
we use the shorthand notation (J3, a) = (fi, a\, . . . , oc m -\) e GF(2) x GF(2) m ~ [ to denote the 
eigenvalues of stabilizer generators. Here e GF(2) is the eigenvalue of the operator S o, namely, 
the phase error detected; and a t e GF{2) is the eigenvalue of the operator Su namely, the bit flip 
error detected, for 1 < i < m - 1. We sometimes abuse the notation to denote a state by (J3,a). 
That is, we denote the states (|0 8m > + |l® m »/ V2 and (|0 8m > - |l® m »/ V2 by (J3,a) = (0,0) and 
(J3, a) = (1,0) respectively. 



B. Depolarization to the GHZ-basis diagonal states 

The m players can depolarize each copy of their shared GHZ state into the GHZ diagonal basis 
using local operation and classical communication (LOCC) in the following way yfl. A player 
randomly chooses an operator from the span of the set of stabilizer generators of the GHZ state 
and broadcast his/her choice to the other players. Then they collectively apply the chosen operator 
to the GHZ state. Since all stabilizer generators of the GHZ state in Eq. © are tensor products of 
local unitary operators X,- or Z, , the players can apply the operator chosen above to the state locally 
using LOCC. Then, they forget which operator they have chosen. The resultant state is diagonal 
in the GHZ basis. Moreover, the error rate of the GHZ state is unchanged in this process. So, we 
can always assume that each state shared among the players are diagonal in the GHZ basis. 

Among all GHZ-basis diagonal states with a fixed error rate (and hence also among all states 
with a fixed error rate), Werner state is the most difficult to work on as far as distillation of GHZ 
states is concerned. This is because one can turn any state into a Werner state with the same error 
rate via a depolarizing channel. Hence, to study the worst case performance of the distillation of 
GHZ states, we suffices to investigate the case in which the input states are Werner states. 



C. Maneva and Smolin's multi-party hashing protocol and its generalization by Chen and Lo 



Maneva and Smolin [4] proposed a multi-party hashing protocol by generalizing the bilateral 
quantum XOR (BXOR) operation |5|] to the multipartite case. In their scheme, the m players 
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carefully choose two (classical) random hashing codes, one to correct spin flip errors and the other 
to correct phase errors, and apply them to their shared noisy GHZ states. This can be done by 
using local operation plus classical communication with the help of a few multi-lateral quantum 
XOR (MXOR) operations. Recall that MXOR is a linear map transforming the state (^)™ \ju h)i 
to ^"j = Q \ji,ji + kj)i for all _/,-,£,• e GF(2). Here, quantum particles with subscript i belong to the 
z'th player. Suppose the source and target states are eigenstates of the stabilizer generator in Eq. © 
with eigenvalues (fii,a{) and (/J 2 , cm) respectively. Then after the MXOR operation, the resultant 
state is also an eigenstate of the stabilizer generator with eigenvalues Q4£| 



MXOR[08 1 ,af 1 ),08 2 ,ar 2 )] = [(ft + j8 2 ,ari),G8 2 ,ar 1 - or 2 )] . (5) 

Using the observation that spin flip error occurred in different qubit of a GHZ state can be 
detected and corrected in parallel, Maneva and Smolin showed that the asymptotic yield of their 
hashing protocol in the limit of large number of shared noisy GHZ states is given by D4J] 

A = 1 - max [H(bd] - H(bo) , (6) 

\<i<m-l 

where H(x) = - £ ,• Pj l°g2 Pj ls m e classical Shannon entropy function. Here the n-bit string bo 
represents the random choice of fi\, . . . where p c corresponds to the eigenvalue of the operator 
So of the Mi GHZ-state |<D m+ )'s and the A^-bit string b { represents the random choice of a u , • ■ • , a m 
where a ei is the eigenvalue of the operator S , of the iih. GHZ state for 1 < i < m — 1 . That is to 
say, H(b ) is the averaged phase error rate and H{bi) is the averaged bit flip rate corresponding to 
the stabilizer generator 5, (for i = 1,2, . . . ,m - 1) over the Af GHZ states respectively. 

Recently, Chen and Lo improved the above random hashing protocol by exploiting the cor- 
relation between the string b\. Specifically, they replaced the spin flip error-correction random 
hashing code used in Maneva and Smolin protocol by the following scheme. (Actually, they only 
considered the case of three players. What we report below is a straight-forward generalization to 
the case of m players as we need to use this generalization later on.) Player 1 applies (classical) 
random hashing to correct spin flip error occurred in his/her share of the GHZ states. He/She then 
broadcasts his/her hashing code used and his measurement results. For i = 2, . . . , m - 1, the z'th 
player carefully picks his/her (classical) random hashing code to correct spin flip error occurred in 
his/her share of the GHZ states based on the broadcast information of players 1 , 2, . . . , i - 1 . Then, 
the z'th player broadcasts his/her code used and his/her measurement results. In this way, the yield 



of Maneva and Smolin scheme can be increased to 
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D 2 = l-max{H(bi),H(b 2 \bi),H(b 3 \b 2 ,bi),--- ,H(b m -i\b m - 2 ,b m - 3 , . . . M)} 

-H(b ) + I(b Q ;b bi) (7) 

where the function /( ; ) is the mutual information between the two classical random variables 
appear in its arguments. 

Applying the random hashing method of Maneva and Smolin to a collection of identical tripar- 
tite (that is, m = 3) Werner states in Eq. ©, one can obtain perfect GHZ state with non-zero yield 
whenever the fidelity F > 0.8075 [|4|]. Using the Chen and Lo's formula in Eq. ©, one can push 
this threshold fidelity down to 0.7554 Q6Q . 



D. Shor-Smolin concatenation procedure and its generalization to the multipartite situation 

n 

Built on an earlier work by Shor and Smolin [1], DiVincenzo et al. introduced a highly error- 
tolerant way of distilling shared Bell states by means of a concatenation procedure |2D. This 
procedure can be generalized to distill shared GHZ states in a straight-forward manner. We report 
this generalization below since we have to use a few related equations later on. 

Suppose m players share Nn copies of imperfect GHZ states for iV » 1 . They perform the 
following two level decoding procedure. First, the players randomly divide these GHZ states 
into N equal parts. Then each player applies a decoding transformation associated with an addi- 
tive [n, k\, d\\ code to his/her own qubits followed by the error syndrome measurements. Surely, 
this can be done with the help of a few MXOR operations. By comparing the difference in 
player's measurement results, they obtain the syndrome 3 e GF(2) {m ~ l ^ n ~ kl) . To continue, each 
party applies another decoding transformation corresponding to a (classical) random hashing code 
[Nki,k 2 ,d 2 ] to correct errors in the GHZ diagonal basis and broadcast the measurement results. 
Finally, they apply the necessary unitary transformation according to the measured error syndrome 
of this random hashing code to get the purified GHZ states. 

Suppose that an additive code [n, ki,d\] is applied and the remaining states after the decoding 
transformation and measurements are denoted by (6, 7) = TRAN [(J3\ , a 1 ) , (J3 2 , a 2 ) , . . . , (fi kx , ou Jl . 



Then, the yield of this concatenated scheme is given by the so-called Shor-Smolin capacity jl 



2] 



D ss = -(l-S x ), (8) 
n 
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where 

s x = J] Pr ^ Kims, rM ■ (4 r) e gf(2)"»}) (9) 

?eG/ ;- (2)(" , -" ( "-*i' 

is the average of the von Neumann entropies of the quantum states conditional on the measurement 
outcomes. Note that in Eq. ©, Pr(s) is the probability that the measurement outcome is s, 

h({pi}) = -J^Pilog 2Pi , (10) 

i 

and 

J> = 1- (ID 

i 

By applying the above procedure to depolarized Bell states using a 5-qubit cat code as the 
inner code and a random hashing code as the outer code (that is, the case of m = 2 and n = 5), 
DiVincenzo et al. found that one can attain a non-zero capacity whenever the channel fidelity 
F > 0.80944 B20. Since the performance of this scheme exceeds that of quantum random hashing 
code and that the 5-qubit cat code is degenerate, the power of using degenerate quantum code in 
quantum information processing is demonstrated. 



E. Other hashing and breeding schemes 



Several other multipartite hashing schemes have been studied |4j,|7|]. In particular, Maneva and 
Smolin's hashing scheme can distill shared GHZ states from copies of Werner states with fidelity 
F > 0.7798 in the limit of arbitrarily large number of players (that is, when m — » oo) M\. An- 
other approach is to use the so-called stabilizer breeding. In particular, Hostens et al. showed that 
stabilizer breeding is able to purify depolarized 5-qubit ring state with fidelity F > 0.756 080. A 
few authors also studied the distillation of graph state subjected to local Z-noise D90 and bicol- 
orable graph state. Jioll Furthermore, Glancy et al. generalized the hashing protocol of Maneva 
and Smolin to purify a much larger class of output state. [1 1] The second column in Table|I] sum- 
marizes the state-of-the-art one-way purification schemes to distill depolarized GHZ states before 
our work. 
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m 


prior art 


our best protocol 


lower bound 


2 


0.8094 


0.8094 


0.7500 


3 


0.7554 


0.7074 


0.6111 


4 


0.7917 


0.6601 


0.5500 



TABLE I: The threshold fidelity of the depolarizing channel above which a GHZ state can be distilled 
by prior art and by our protocols. Also listed is the lower bound of the fidelity below which no one-way 



protocol can distill shared GHZ state using Eq. (1431 in Sec. IIII Ei As for prior art, the threshold fidelity for 
m = 2 is given by the 5-qubit cat code |2D. For m = 3 case, the threshold is computed by the Chen and Lo's 
formula in Eq. ([7]). For m = 4, the threshold is given by the Maneva and Smolin's hashing protocol (4J] 
in Eq. ©. 



III. OUR PROTOCOLS INVOLVING DEGENERATE CODE AND THEIR PERFORMANCES 



A. Our protocols 

Our three protocols are natural extensions of the Shor-Smolin concatenation procedure to the 
case of purifying GHZ states. They all use the same degenerate quantum code as the inner code. 
Specifically, suppose the m players share Nn copies of Werner state with N » 1 . As shown in 
Fig. [TJ to distill perfect GHZ state, each player applies the (classical) [n, l,n] repetition code, 
whose stabilizer generators are 

ZoZi , Z0Z2 , . . . , ZoZ„_ 1 , (12) 

to his/her own n qubits. That is to say, they randomly partition the Nn shared noisy GHZ states into 
N sets, each containing n noisy GHZ states. In each set S, they randomly assign one of the noisy 
GHZ state as the source (and call it the 0th copy of |O m+ ) in the set) and the remaining (n - 1) noisy 
GHZ states as the targets (and call them the j'th copy of |O m+ ) in the set for j = 1,2, . . . ,n - 1). 
They apply the MXOR operation to copies of |O m+ ) in each set and then measure all the target GHZ 
states in the standard computational basis while leaving all the source GHZ states un-measured. 
We denote the syndrome and the remaining state in each set by s* s e Qp(2f n ~ l ^ m ~ l) and (6 3 , y s ) 6 
GF(2) m respectively. (Since the partition into sets is arbitrarily chosen and our subsequent 
analysis only makes use of the statistical properties of the states in each set, we drop the set label 
S in all quantities to be analyzed from now on.) 

Our three protocols differ in the use of outer codes. For the first protocol, each player applies 
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FIG. 1: Our three GHZ state distillation protocols for m = n = 3. They all use the [n, l,n] repetition code 
as their inner codes; and they differ by the kind of random hashing outer code used. 

a (classical) random hashing code [N,k 2 ,d 2 ] that corrects GHZ diagonal basis errors to the N 
remaining states (each coming from a different set S) and exchanges the measurement results. 
Clearly, this protocol is reduced to the Shor-Smolin concatenation procedure [1] when m = 2. 

For the second protocol, the players follow Maneva and Smolin's idea i§\ by using two (classi- 
cal) random hashing codes, one to correct spin flip error and the other to correct phase shift. In this 
sense, the outer code used in our second protocol is a random asymmetric Calderbank-Shor-Steane 

n 

(CSS) code. Whereas players in our third protocol use the Chen and Lo's generalization [6] as their 
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outer code. That is, the outer code is a random asymmetric CSS code whose decoding circuit is 
carefully designed to exploit the correlation between the bit string b t . 

In all the above three protocols, the players have to apply the corresponding unitary transfor- 
mation for the outer code to obtain the purified GHZ states. (See Fig. CD) 

Clearly, the yield of the first protocol is the Shor-Smolin capacity given by Eqs. © and ©. 
And by applying Eq. © to the noisy GHZ state to be fed into the outer code of the second protocol, 
we conclude that the yield of the second protocol equals 



D 



1 



MS 



1 - 

V ?eGf(2)('"- i;! ' - h 



V Prtf) ( max [H( 7i \t)\ + H(6\^) 

W <i<m- 1 ) 



where 



H<yi\f) = h({ Pr « 6 > ?M : * e GF ™ 

S,yi ,...,7i-i ,Ti+ 1 •■ ■ ..7m- 1 sGF(2) 

for i = 1,2, . . . ,m - 1, and 

H(m = h({ J] Pr((<5,r)|^):<5eGF(2)}). 

yeGF(2) m - 1 

Similarly, from Eq. ©, the yield of the third protocol is 



(13) 



(14) 



(15) 



1 



D C L = - 



1 



V Pr(^) ( max [Hiy^uji-i, ■ . • , 7u ft] + H(S\2) - 1(6; Y\?>) (16) 

^— ' (\<i<m-l ) 

V seGF{2f m - l " n - l '> J 

where 

H(y i \y i - l ,yi-u...,yi,f) = K{ J] ms.r)fyi-u7i-2,...,yi,3) ■ Yi e GF(2)}) (17) 

S,y i+ i,y i+ 2,-;y m -teGF(2) 

for i = 1, 2, . . . ,m — 1 and 7(5; yls) is the conditional mutual information between 6, and 7 given s. 



B. Evaluating the yields for Werner states for our protocols 

To analyze the performance of our three protocols when applied to Werner states, we first have 
to calculate the distribution of outcomes after passing the Werner states through the inner repetition 
code. For an arbitrary but fixed set S, using the compact notation introduced in SecUH we denote 
the error experienced by the j'th copy of |O m+ ) in this set by (fij, aj) for j = 0, 1, . . . , n - 1. After 
decoding the inner code, namely, the [n, l,n] repetition code whose generators of the stabilizer are 
written down in Eq. (fT2l) . the syndrome = s 2 , ■ ■ ■ s n -\) 6 GF(2) (m " 1)( " _1) obtained obeys 

sj = aj - a e GF(2) m - 1 (18) 
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for all 1 < j < n — 1. Furthermore, the remaining state shared among the m players is 

n-l 

(S,y) = (Y J Pj,a ). (19) 

7=0 

To simplify notation in our subsequent discussions, we define 

so = (20) 

so that Eq. (PT8T) is also valid for j = 0. 

To evaluate the capacity for each of our three protocols, we first have to calculate the conditional 
probabilities Pr((£, y)\it) and Pr((£, y)\y\, y 2 , • • • , Ji-i, st) in Eqs. © and (fT6l) . respectively. We begin 
by computing the probability Pr((£, 7) Al) that the source state has experienced the error (5, y) = 
(Z"=o/3> Q'o) after the decoding transformation of the repetition code in Eq. (Q3 and that the error 
syndrome for the repetition code is s € GF(2) (m " 1)( "" 1) . Clearly, 

Pr((5,y)A^ = Pr((£(sUy)) (21) 

where 

n-l 

<£(Z 8, y) = {(fij, aj) n rl e GF(2) mn -.5 = ^^,? = a , a e = s e +a for t = 1, 2, . . . , n-l} . (22) 

7=0 

Since the repetition code and our decoding transformation are highly symmetric in the sense 
that they are invariant under relabeling of qubits, it is not surprising that the set (£(s*, 6, y) is in- 
variant under permutation of phase errors. That is to say, (J3j, of ; )"~o 6 ®C^>£y) if and only if 
(J3 n (f), a j)"Zo 6 $ 7) where ^ is a permutation of {0, 1, . . . ,n - 1}. 

1. Finding Pr{{8, y) A f) 

We proceed by introducing the concepts of depolarization weight and depolarization weight 
enumerator similar to the ones proposed by DiVincenzo et al. D2j]. Let (J3j, aj) be the state of the 
jth noisy GHZ state shared among the m players. The depolarization weight of the order n-tuple 
(fij, a j)"Zo 6 GF(2) mn is defined as its Hamming weight by regarding this rc-tuple as a vector of 
elements in GF(2)" 1 . In other words, 

wt ((fij, aj) n f^) = \{j G {0, 1, . . . ,n - 1} : (fij, aj) * (0, 0)}| . (23) 
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Physically, the depolarization weight measures the number of shared GHZ states that experienced 
an error; thus, it is invariant under permutation of the n possibly imperfect GHZ states. Since a 
GHZ state has equal probability of having each type of error after passing through a depolarizing 
channel, there is an equal probability for the n depolarized GHZ states to experience errors with 
the same depolarization weight. Thus, we may find the probability Pr((5, 7) A s) by studying the 
depolarization weight enumerator w( (£(,?, 6, 7); x, y) where 

w{%; x,y) = Yj x wt( V" wt(3) • (24) 

The depolarization weight enumerator of a set is a natural generalization of the concept of weight 
enumerator of a code. 

Finding an explicit expression for the above depolarization weight enumerator for an arbitrary 
set or coset is a very difficult task. It is the high degree of symmetry in the repetition code that 
makes this task possible. In fact, one may transform a state in 8, 7) to another state in the same 
set by applying phase shifts to a few qubits. 

By counting the number of different possible combinations of (J3j,aj)'s subjected to the con- 
straint that (J3j, (XjTjll e (£(s, 6, 7), we have 

1! 

.1 



wWZ 8, 7); x, y) = Y' -j=^ (25) 

[ I V 



>,t- 

>£GF(2), 
(eGF(2)™- 1 

where the primed sum is over all a t /s satisfying the constraints 

a u t > Vi, t , (26) 
J] a u = n , (27) 

ieGF(2), 
(eGF(2) m -! 

2 a u = \{je{0,h...,n-l}:s j + 7 = t}\\/t (28) 

;eGF(2) 

and 



J] i<*i,t = 8. (29) 



ieGF(2), 
teGF(2) m - { 

Note that the symbol a tJ in the above equations can be interpreted as the number of GHZ states 
that experienced the error (z, t) before the commencement of our distillation protocol. 
Let 

k = k(l 7) = \{j g {0, 1, . . . , n - 1} : sj + 7 * 0}| (30) 
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be the number of qubits having spin flip for each element in (£(s, 6, y). We have two cases to 
consider. 

Case (a) k > 0: That is, there exists i such that S[ + y ± 0. Hence, wt((/? ; -, Sj + y)pj) is independent 
of the value of /3 e e GF(2). In addition, by regarding the equation Z"=oA/ = 5 as a. bijection 
relating fi e e GF{2) and 6 e GF(2), we conclude that the depolarization weight enumerator 
w(<&(&, 6, y); x, y) is independent of the value of 6 e GF{2). Hence, 

wmi6,y);x,y) = ]- Y" (31) 

n 

(eG/ r (2) m - 1 

where the double primed sum is over all a ( /s satisfying constraints Eq. (|26l)-(|28l) only. Conse- 
quently, 

wW,6,y);x,y) = \ Y — f "*) x^V ' 

2 fej n a ^ 1 aa ° ; 



(eGF(25 m-1 \(0) 



^«-«o,<yo,o 



(n — k\ 

, a o,o 

«0,0 

= 2*- 1 x*(x + y)"-* . (32) 



Case (b) = 0: That is, Sj ■ + y = for all j so that phase shift is the only type of error a GHZ state 
may experience. In this case, the union of disjoint sets U<5eGF(2) 7) i s equal to the set of all 
possible phase errors experienced by the n shared GHZ states. As a result, 

£ wWZ S, 7); x, y) = w( [j m 6, y)-x,y) = Y i U *Y = (* + 7)" • (33) 

<$eGF(2) SeGF(2) i > Z ' 

Similarly, 

2 (-1^(^,7);^) = ^ |J ^47);-^) = ^ " (-x)y-'' = (y-xr. (34) 

<SeGF(2) <SeGF(2) i 

(Note that the validity of the above equation follows from the observation that 6 = 1 if and only 
if the number of qubits having phase shift error before the commencement of our protocol is odd. 
Moreover, the coefficient of w({a}; -x,y) is negative if and only if wt(a) is odd.) From Eqs. (1321)- 
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(1341) . we conclude that 



w(&(f, 6, 7); y) = w(k, 6;x,y) = 



2 k - l x k (x + y) n - k 



if < k < n, 



- [(x + yf + (y - x) n ] if k = and 8 = 0, (35) 



2 K* + ~ (y ~ *)"] if ^ = and 5 ^ 0, 

where k = k(s, 7) is the number of GHZ states that experienced some kind of spin flip for each of 
the state in (£(s, 6, 7) as defined by Eq. (|30l) . 

Recall that (£(?, 6, 7) is invariant under permutation of phase errors among the n GHZ states. 
Moreover, both the depolarization weight and the value of k{^,y) are invariant under permu- 
tation of the n GHZ states. So, it is not surprising that the depolarizing weight enumerator 
w((£(s, 6, 7); x, y) depends only on the values of k and 6. Therefore, our shorthand notation 
w(k, 8; x) makes sense. 

From Eq. (ED and by substituting x = (1 - F)/(2 m - 1), y = F into Eq. ([35]), we find that 

( i k ~\\ - F) k (2 m F -2F+ l)"- k . p n , 

if < k < n, 



Pr((<S,7)AsO 





(2 m - 


1)" 




(2 m F - 


2F + 1)" 


+ (2 m F - 


1)" 




2(2 m - 


1)" 




(2 m F - 


2F + 1)" 


- (2 m F - 


1)" 




2(2"' - 


If 





if k = and 5 = 0, 



if jfc = and 5^0. 



(36) 



for a depolarizing channel with fidelity F. Note that by fixing the number of players m, the number 
of noisy GHZ states shared between the players n and the fidelity of the depolarizing channel F, 
the probability Pr((5, 7) A it) can take on at most (n + 2) different values. 



2. Finding Pr((6, y)\s) 



Clearly 



Pr(s) = Pr(sAf) 

teGF(2f'- 1 



(37) 
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where Pr(s At) is the probability that the error experienced by the n noisy GHZ states is (J3j, Sj+t)"^ 
for/3 y e GF(2). For a depolarizing channel with fidelity F, 



teGF(2f 



2(1 - F) 
2 m - 1 



k(s,t) 



F + 



l-F 
2 m - 1 



n-k(s,t) 



(2 m - 1) 



^— £ /K02«"(l - F)'(2 m F -2F + 1)"-' 



i=0 



where 



Therefore, 



/K0 = |{i e GF(2r _1 :^,0 = i}|. 



n—k 



2*-i(l - Ff(2 m F - IF + 1) 
Z,/?0')2<(1 - F)''(2 m F - 2F + l)"- ; 

(2»F - 2F + l) n + (2 m F - l) n 
2 Zi/K02'(l - F)''(2 m F - 2F + 



(2 m F-2F + l)"-(2 m F- 1)" 



if <k<n, 



if it = and 5 = 0, 



if it = and 5 ± 0. 



(38) 



(39) 



(40) 



I 2 Z,-/K02'(l - F)'(2 m F -2F+ I)"-' 

So combined with Eqs. © and ©, we have a working expression for /z({Pr((5, y)|s) : (5,7) € 
Gi 7 (2) m }) and hence Z) S s- While combined with Eqs. (fT3l-(fT5l). we have a working expression 
for h{yi\i), h{5\?) and hence D MS . In the calculation of A(yil7i-i>7i-2> •■•,7i ! s), we need to first 
compute the probability Pr((5, y)l7,--i, 7,-2, • • • ,71, As for repetition code, Eq. (PT8l tells us that 
the kind of spin flip error experienced by the jth GHZ state aj is known once 7 and s are fixed. 
Hence, Pr((5, 7)|y,-i , 7,-2, . . . , 71 , s) is also given by Eq. (l40l) . Consequently, using Eqs. (fT5l)-([T7l). 
we get a working expression for h{yi\yi-\, ...,y\,f), h(5\s), 1(6; y\s) and hence D C l- 



3. Complexity issue on the computation ofD$s, Dms an d Dcl 

Apparently computing D ss , D MS and D CL using Eqs. ©, @, (fT3l)-([T71). (I3~8l)-(l40l) are ex- 
tremely inefficient as the sum on ^may take on 2 m(,I_1) possible values. Nonetheless, the numerical 
values of many terms in the R.H.S. of Eq. © are the same because the if dependence of Pr(s) and 
Pr((5, y)\s) come indirectly from the distribution of {k(s, t) : t e GF{2) m ~ 1 }. Note that there are at 
most Yn =<i Pi(n) different possible distributions for {k(s, t) ; t G GF(2) m ~ 1 } where Piin) denotes 
the number of ways to express n as a sum of exactly i positive integers. Moreover, P\(n) scales 
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as exp(7r yJ2n/3)/4n V3 in the large n limit II 1211 . Consequently, for a fixed m, we may regroup the 
sum Eq. © so as to compute S x by summing only sub-exponential in n terms. Although this is 
not a polynomial time in n algorithm, it is good enough to obtain the numerical values for S x and 
hence the yield of our first protocol, namely, the Shor-Smolin capacity D ss for a reasonably large 
number of n. By the same token, the yields of our second and third protocols, namely D MS and 
D C l respectively, can also be computed in sub-exponential time in n. 



C. Performance of our three schemes 



We study the performance of our three protocols by studying the yield as a function of the 
channel fidelity F . In particular, we would like to find the threshold fidelity, namely, the minimum 
fidelity above which D > 0, as a function of the number of players m and the repetition codeword 
size n. And we denote the threshold fidelities for our first, second and third protocols by F^Jm, n), 
F^ n (m, n) and F ^J jn, n), respectively. 



1. Subtlety in the computation of threshold fidelities 

Finding the values of F ^J m, n), F M . s (m, n) and F^ n (m, n) requires extra care. Let us explain 
why for the case of F^ n (m, n). And the reason for the other two cases are similar. 

Since S x is a continuous function of the channel fidelity F, Eq. © implies that F^ m {m,n) is 
the root of the equation Sx = 1- Note that 

1 -S x = Pr(0)[l -h({Pr((6,y)\6) : (8, y) e GF(2) m })] 

- ^ Pr(s)[/*({Pr((<S, yM : (6, y) e GF(2) m }) - 1] . (41) 

From Eq. (00]), we know that for F » 1/2, /z({Pr((£, y)\i) : (6, y) e GF(2) m }) is less (greater) than 
1 if s = (s + 0). More importantly, for a fixed m, lim,,^ h({Pr((6,y)\s) : (6, y) 6 GF(2) m }) = 
1~(1 + ) for s = (s ± 0). Thus, Eq. (|4TT) shows that 1 - S x is the difference between two small 
positive terms. This makes the computation of F^ n (m, n) together with the analysis of its trend 
as a function of m and n, particularly for a large n, difficult. Even worse, for F < 1 and for a 
sufficiently large n, the errors experienced by the noisy GHZ states (6, y)" =() satisfying s = are 
not in the typical set. Actually, we found that for F close to F^ n (m, n), the dominant terms in the 
R.H.S. of Eq. (|4TT) almost always correspond to atypical errors experienced by the GHZ states. 



16 



In spite of these difficulties, we are able to accurately compute the yield of our first protocol 
£>ss, namely, the Shor-Smolin capacity, as a function of the channel fidelity F for the classical 
[n, l,n] repetition code acting on the |<D m+ )'s. And from this, we can deduce the correct threshold 
fidelity for our first protocol F^Jm, n) as a function of the number of players m and the number of 
shared noisy GHZ states n. The trick is to use rational number arithmetic to obtain an expression 
for S x for a given rational number F before converting this expression to an approximate real 
number. The same trick also enables us to obtain accurate values for F ™? (m, n) and F^rjm, n), 
namely, the threshold fidelities of our second and third protocols. 



2. The superior performances of our three protocols 

The yields of our three protocols are shown in Figs. [2HS and the corresponding threshold 
fidelities are tabulated in Tables HTj-HVl By comparing these tables with the second column of 
Table HI we make the most important conclusion of this paper: for the multipartite case (m > 3) 
and for any number of shared GHZ states n, the error-tolerant capability of our third protocol is 
strictly better than our second, which is in turn strictly better than that of our first. And under the 
same conditions, the error-tolerant capability of our first protocol is already better than the best 
scheme in literature before this work. So once again, we show the powerfulness and usefulness of 
degenerate codes in one-way entanglement distillation. 

Whereas for the bipartite case (m = 2), Tables HH-HVl show that all our three protocols can 
tolerate almost the same level of error. It means that the use of random asymmetric CSS outer code 
does not give any significant advantage here. (Actually, we find that using random asymmetric 
CSS outer code decreases the error-tolerant capability for n < 3.) Interestingly, the threshold 
fidelities for our second and third protocols agree to at least four significant figures for m = 2. 
This finding can be understood as follows. As we have discussed in Sec. IIII Bl and particularly in 
Eq. (1401) . the probability of 5 = equals the probability of 8 = 1 provided that s + irrespective 
of the value of y. That is to say, 1(6; y|S) = whenever 0. So, it is not surprising to find that 
the weighted mutual information }Zs Pr(^)/(5; y\ft) becomes negligibly small when the fidelity F 
is close to its threshold value. Combined with Eqs. (fT3l) and (fT6l) . it is not unnatural to find that 
F MS (2,n) * F CL (2,n). 

min v ' ' min v ' ' 

Our numerical computation shows that lim^o, F^Jm, n), lim,,^ F^fj jn, n) and 
lim„^oo F CL (m, n) are decreasing functions of m. Besides, lim^oo F ss (4, n) is smaller than 
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0.8098 


0.7582 


0.7185 


0.6898 


0.6696 




11 


0.8104 


0.7492 


0.7021 


0.6677 


0.6435 




15 


0.8110 


0.7449 


0.6938 


0.6565 


0.6301 




21 


0.8118 


0.7416 


0.6870 


0.6471 


0.6188 




31 


0.8128 


0.7391 


0.6814 


0.6390 


0.6089 


TABLE II: The threshold fidelity F ss as a function of m and n. 
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m 
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3 


4 


5 
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2 


0.8137 


0.7788 


0.7541 


0.7369 


0.7253 




3 


0.8101 


0.7631 


0.7261 


0.6991 


0.6781 




4 


0.8102 


0.7551 


0.7091 


0.6781 


0.6571 


n 


5 


0.8095 


0.7566 


0.7111 


0.6771 


0.6521 




6 


0.8100 


0.7522 


0.7081 


0.6721 


0.6421 




7 


0.8098 


0.7501 


0.7051 


0.6711 


0.6441 




11 


0.8104 


0.7475 


0.6951 


0.6581 


0.6311 




15 


0.8110 


0.7446 


0.6901 


0.6511 


0.6221 



TABLE III: The threshold fidelity F MS as a function of m and n. 



0.7798, the fidelity threshold of the Maneva and Smolin's hashing scheme in the large m limit [4J]. 
So, for m > 3, our three protocols all tolerate a higher noise level than all other one-way schemes 
known to date. Figs. [2H1] further depict that the yields of our protocols D ss , D M s and D C l are very 
steep functions of F around their corresponding threshold fidelities. Thus, a reasonable yield can 
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n 
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0.8095 
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0.6776 


0.6582 
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0.8100 


0.7099 


0.6684 


0.6551 


0.6217 
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0.8098 


0.7086 


0.6650 


0.6480 


0.6133 




11 


0.8104 


0.7081 


0.6642 


0.6372 


0.6062 




15 


0.8110 


0.7074 


0.6601 


0.6284 


0.6036 



TABLE IV: The threshold fidelity F Ch as a function of m and n. 

J mm 



be obtained when F is equal to, say, 0.02 higher than the threshold. 

Another interesting feature found in Tables ITTl-llVl and Figs. [2H3]is that the threshold fidelities 
F ss (m,n), F MS (m,n) and F CL (m,n) are all decreasing functions of n for m > 3. That is, our 

min v ' 7 ' min v ' ' min v ' ' ° ' 

protocols attain a higher capacity if players use a longer repetition code whenever m > 3. In 
contrast, DiVincenzo et al. found that F^J2,n) attains global minimum when n = 5. Besides, 
F ^ n (2, n) > F ^ n (2, n+l) for a small even integer n [2J]. Interestingly, Tables Hill and ITVl show that 
F ^f n {2, n) and F^ n (2, n) behave in the same way, too. 

Lastly, we remark that for m > 3, the improvement in the error-tolerant capability for increasing 
n comes with a price. For a fixed m > 3, the yields of our protocols decrease as n increases 
provided that the channel fidelity F is close to 1 as depicted in Figs.[2H4l This is because as n 
increases, more shared GHZ states must be wasted in order to obtain the error syndrome s even if 
the channel is noiseless. 



D. Understanding the trend of the threshold fidelities 

Although the discussions in this subsection focuses on the trend of the threshold fidelity of our 
first protocol, namely, F^ n (m, n), the essential ideas also apply to the cases of our second and third 
protocols, that is, F™?(m, n) and F^ n (m, n). 

The reason why F^ n (2, n) is a sawtooth- shaped function of n for n < 8 is related to the behavior 
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FIG. 2: The yield Dss of our first protocol for distilling |cD" ,+ ) after passing through a depolarizing channel 
of fidelity F using the classical repetition code [n, l,n] as the inner code for various n when (a) m = 2, 
(b) m = 3, (c) m = 4 and (d) m = 5. 

of h({Pr((6,yM : (6, y) e GF(2) m }). It is easy to check that for m = 2, fc({Pr((cS, : (5,7) e 
Gi 7 (2) m }) is equal to (much less than) 2 provided that the depolarization weight wt(s) = n/2 
(wt(s) ^ n/2). For a small even n, there is a non-negligible probability of finding s with wt(s) = 
n/2 so that the root of S x = 1 and hence the value of F^ in are determined mainly by the summing 



20 




FIG. 3: The yield Dms of our second protocol using the same sets of parameters as in Fig. [2] 

only over those s's with depolarization weight or n/2 in Eq. (|4TT) . In contrast, for a small odd 
n, all entropies in the R.H.S. of Eq. (|4TT) are much less than 2. Hence, the corresponding value of 
F^J2, n) is lower than F^J2, n + 1). In other words, the reason for F^Jl, n) > F^J2, n + l) for 
a small even n is that there is a non-negligible chance that exactly half of Bell states used by the 
inner repetition code have spin flip error so that players have absolutely no idea what kind of error 
the remaining unmeasured Bell state has experienced. 
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FIG. 4: The yield Dql of our third protocol using the same sets of parameters as in Fig. [2] 

However, the situation is very different when m > 3. In this case, the condition for 
/?({Pr(((5, y)\s) : (6, y) e GF(2)' n }) > 2 is that one can find an integer i such that fg(i) > 2 and 
fs(j) = f° r all j < i. More importantly, for a depolarizing channel with F > 1/2, the probability 
Fr(s) of finding this kind of s with /z({Pr((£, y)\s) : (6, y) e GF(2) m }) > 2 is much less than that in 
the situation of m = 2. Thus, the contribution of terms with entropy greater than or equal to 2 in 
Eq. (HTI) becomes much less significant when m > 3. So, it is not surprising to find that for a fixed 
in > 3, F ^. n (m, n) is not a sawtooth-shaped function of n when n is small. 
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It is also easy to understand why lim,,^ F^(m, n) is a decreasing function of m. One simply 
check by Taylor's series expansion that in the limit of large n and for a fixed 1/2 < F < 1. Then 
we find that the first term in the R.H.S. of Eq. (|4Tj) is an increasing function of m; and that the 
summand in the second term in the R.H.S. of Eq. (|4T|) is almost surely a decreasing function of m 
in the large n limit. 

As we have pointed out that the value of F ^ n (m, n) depends on the entropy of a few atypical 
set of errors experienced by the GHZ states. We do not have a good explanation why F ^ n (m, n) is 
a decreasing function of n for m > 3. 



E. Breaking the F > 0.75 limit? 



No t error correcting quantum code of codeword size At exists |5ill3j]. Hence, it is impossible 
to distill Bell states using an one-way scheme provided that the fidelity of the depolarizing channel 
is less than or equal to 0.75 150 . That is why F^(2, n), F^(2, n), F^(2, n) > 0.75. Interestingly, 
a few F ss (m, n)'s, F MS (m, nYs and F CL (m, n)'s listed in Tables iHWTVl are less than 0.75. Does it 

min v 5 7 ' min v ' 7 mm v ' 7 1 — 1 1 1 

make sense? 

To solve this paradox, let us recall that the Pauli errors experienced by a GHZ state shared 
among m players can always be regarded as taken place in (m - 1) of the m qubits. From Eq. ©, we 
may regard that at most one of the (m-1) qubits may experience a phase error. So, the probability 
that a depolarized GHZ state has experienced phase error but not spin flip is ( 1 - F) /(2 m - 1 ), where 
F is the channel fidelity. And the number of erroneous qubits equals 1 in this case. Besides, the 
probability that exactly i out of the (m-1) qubits have experienced spin flip is 2(1 -F)( m T 1 J/(2 m - 1) 
for i = 1, 2, . . . , m - 1, where the extra factor of 2 comes from the fact that the spin-flipped GHZ 
state may experience phase shift as well. Hence, the average number of erroneous qubits divided 
by (m - 1) is given by 



1 



e 



1 1 ! 



1 



m-1 



1 - F 2(1 - F) . m-1 



1 



+ 



1 - F 

2 m - 1 



>m— 1 



+ 



1 



m 



1 



(42) 



130, g< 1/4. 



Since no t error correcting quantum code has codeword size less than or equal to At |5|, . 
Consequently, a lower bound for the threshold fidelities F* in (m, n) (for x = SS, MS, CL) is given 
by 



Fl-(m,n) > Fbomid = 1 - 



1 



\tn-\ 



1 



+ 



m 



1 



(43) 
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A quick look at the third and the fourth columns in Table U convinces us that our protocols do not 
violate this general limit. Actually, one of the reasons why we can distill shared GHZ states when 
F < 0.75 for m > 3 is that the average qubit error rate for a depolarized GHZ state is given by 
Eq. ([42]) . which is smaller than (1 - F). Note in particular that in the large m limit, the average 
qubit error rate for a depolarized GHZ state is close to 1 /2. So, it is not surprising that the bound 
Abound approaches 1 /2 in this case. 



IV. GENERALIZATION TO HIGHER DIMENSIONAL SPIN 



A. Our extended protocols 

Our three protocols can be generalized to the case when the Hilbert space dimension of each 
quantum particle is greater than 2. That is to say, the m players wanted to share the state 

'i-i 



through a depolarizing channel by one-way entanglement distillation. The quantum codes used in 
the three generalized protocols are extensions of their corresponding binary codes to the g-nary 
ones. In particular, their common inner code becomes classical [n, l,n] q repetition code. 
We have the following two cases to consider. 

1 . For q = p m where p is a prime number, we may impose a finite field structure GF(q) to the 
system by defining 

Xj:\i)^\i + j) (45) 

and 

Zj : |0 i * oj t ;^ \i) (46) 

for all j e GF(q) where co p is a primitive pth root of unity, Tr is the absolute trace and all 
arithmetic are performed in the finite field GF(q). 

2. Alternatively, for any integer q > 2, we may impose a ring structure Z/qZ to the system by 
defining 

Xj:\i)^\i + j) (47) 

and 

Zj-.\i)^a/j\i) (48) 
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for all j eZ/qZ where co q is a primitive qth root of unity and all arithmetic are performed in 
the ring Z/qZ. 

From now on, we use the symbol K to denote either the finite field GF(q) or the ring Z/qZ. 
Similar to the case of q = 2, we use the compact notation (J3, a) = (fi, a\, a2, . . . , a m -\) to denote 
the eigenvalue of the stabilizer generators where /3 e K and a e K m_1 . 

In the qubit case (that is, q = 2), the error syndrome measurement is performed with the aid of 
CNOT gates. In the case of q > 2, this can be done via the operator \i, j) — > \i, i-j) for all i, j e K. 
Suppose the error experienced by the 7th copy of I®™" 1 ") is (J3j, aj) for j = 0, . . . ,n — 1. Then after 
measuring the error syndrome for the classical [n, l,ri] q repetition code, we get s* = (s\,..., s n -i) 
where 

Sj = aj - a (49) 
for 1 < j < n - 1. Furthermore, the remaining state shared among the players becomes 

n-l 

(6, 7 ) = (J]/3j,a ). (50) 

7=0 



B. Finding the capacities of our three generalized protocols 

The analysis in Sec. IIIIBI can be easily generalized to the case of qudits (that is, q > 3). In 
particular, we prove in the Appendix that 

' q k ~ x {\ - F) k (q m F - qF + q - l) n ~ k 



Pr((6, 7 )\s) = 



Mi) [q(l ~ F)] 1 (q"'F - qF + q - 

(q m F - qF + q - 1)" + (q - 1) (q m F - If 
q iZi m ted - F)]> (q m F - qF + q - 1)""' 



(q m F - qF + q - If - (q m F - 1)" 



if k > 0, 



if k = and 6 = 0, 



if k = and 6 * 0, 



(51) 



v q Xi Mi) [qa ~ F)] 1 (q m F - qF + q - 1)"" 
where k and fg(i) axe given by Eq. (1301) and Eq. (IA9I) respectively. 

The yields of our three generalized protocols can be computed using Eqs. ®-([9]) and (fT3l)-(fT71) 
just like the case of q = 2. Nevertheless, there is an important subtlety. Since the players can make 
full use of each of the q possible error syndrome measurement outcomes to distill the generalized 
GHZ state the entropies and conditional entropies in Eqs. ©, (fT4l) . (fT5l) and (TT71) should be 

measured in the unit of dit rather than bit. That is to say, the base of the logarithm used in these 
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0.7441 


0.6927 


0.6611 


0.6443 


0.6370 
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0.6877 


0.6530 


0.6337 


0.6246 




11 
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0.5962 
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0.7449 


0.6700 
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TABLE V: The threshold fidelity F as a function of m and « for q - 3. 

^ mm 1 

entropies should be q instead of 2. And since the dimension of each information carrier q has 
changed, one should not directly compare the yields of the qudit-based protocols with those of the 
standard qubit-based ones. Note further that similar to the original qubit-based protocols, we can 
compute these yields in a time sub-exponential in n. 



C. Performance of our generalized protocols 

Figs. [5]47] depict the yields of our three generalized protocols in the case of q = 3. And Ta- 
bles IVHYlIl list the corresponding threshold fidelities. Clearly, the trend of the threshold fidelities 
of our three generalized protocols are very similar to their corresponding cases of q = 2. In partic- 
ular, for q = 3, the threshold fidelities F ss (m, n), F MS (m, n) and F CL (m, n) are decreasing function 

' 1 ' min v ' min v ' ' mm v ' ' ° 

of n for any fixed integer m > 3; while they reach global minima at n = 7 provided that m = 2. 
These findings are not completely surprising because the arguments we have used to explain the 
trends of the yields and threshold fidelities for our three protocols in the case of q = 2 reported in 
Sec. lIIIDl are also applicable here after minor adjustments. 
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LE VI: The threshold fidelity F MS 
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5 


0.7444 


0.6337 


0.6104 


0.5916 


0.5895 




6 


0.7443 


0.6334 


0.6061 


0.5896 


0.5855 




7 


0.7443 


0.6304 


0.6023 


0.5877 


0.5829 




11 


0.7446 


0.6249 


0.5910 


0.5790 


0.5670 




15 


0.7451 


0.6235 


0.5865 


0.5701 


0.5560 



TABLE VII: The threshold fidelity F CL as a function of m and n for q = 3 

J min -* 



D. Lower bound for the three threshold fidelities 



The proof that no t error correcting quantum code with codeword size At is also applicable to 
qudits II 1411 . We may use this fact to establish a lower bound for the threshold fidelities of our three 
generalized protocols when qudits are used as information carriers. Since the proof is also the 
same as that of the qubit case reported in Sec. IIII El here we only write down the bound without 
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FIG. 5: The yield Dss of our first generalized protocol for distilling ) after passing through a depolar- 
izing channel of fidelity F using the classical repetition code [n, l,n] q as the inner code for various n when 
(a) m = 2, (b) m = 3, (c) m = 4 and (d) m = 5. 

giving the details of the proof: 

F s ±(m, n), FZ(m, n), F^m, n) > F bound = 1 - |— - L^ 1 + — - . (52) 
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FIG. 6: The yield Dms of our second generalized protocol for distilling |®™ + ). The parameters used are the 
same as those in Fig. [5] 



V. SUMMARY AND DISCUSSION 



In summary, we have introduced three one-way GHZ state purification protocols using degen- 
erate codes by extending the works of DiVincenzo et al. on one-way Bell state purification via 
degenerate codes as well as the works of Maneva and Shor |4j] and its generalization by 
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FIG. 7: The yield Dql of our third generalized protocol for distilling |3>^ + ). The parameters used are the 
same as those in Fig. [5] 

Chen and Lo [6] on multipartite entanglement purification using random asymmetric CSS codes. 
Then, we calculate the yields of our three protocols when the inputs are Werner states. The method 
we used to calculate these yields is divided into two steps. The first step is to calculate entropies 
or conditional entropies such as h({Pr((5, y) : (5, y) e GF(2) m }) by means of the so-called de- 
polarization weight enumerator. Actually, the first step can be easily extended to the case of an 
arbitrary stabilizer inner code, an arbitrary un-correlated noise model and an arbitrary stabilizer 



30 



output state. The second step involves the computation of a weighted sum of the entropies or 
conditional entropies obtained in the first step. Nonetheless, for a general stabilizer inner code, a 
general un-correlated noise model and a general output state, this sum may not be practical as it 
involves up to about 2" m number of terms. Fortunately, as the inner code used in our purification 
scheme is the highly symmetrical classical repetition code, we are able to greatly simplify the sum, 
making the computation of the yields in a time which is sub-exponential in n when the GHZ states 
are subjected to depolarization errors. In this way, we can calculate the corresponding thresh- 
old fidelities accurately and reasonably fast. This is quite an accomplishment because finding the 
threshold fidelities involves the accurate determination of the sign of the difference between two 
small positive numbers provided that the number of players m > 3 and the codeword size of the 
inner repetition code n is large. (See, for example, Eq. (EH)-) Just like the Bell state case, we 
discover that the threshold fidelities of our three protocols are better than all known one-way GHZ 
state purification schemes to date. So, once again, the power of using degenerate codes to combat 
quantum errors is demonstrated. 

We also extended our scheme to tackle the case when the information carriers are qudits instead 
of qubits. We find that the performance trend of these generalized schemes are quite similar to 
those of the qubit cases. 

There are a few un-answered questions, however. Here we list some of them. The reason why 
the threshold fidelities F ss F MS and F CL decrease with n for m > 3 is not apparent. And apart 

mm mm mm r 

from the general statement that degenerate codes pack more information than non-degenerate ones 
making them powerful in one-way purification of GHZ states, can we specifically understand why 
using classical repetition code concatenated with a random hashing quantum code is more error- 
tolerant than a few other choices of degenerate codes? Ill5n Along a different line, it is important to 
find out the value of lim,,,-^ lim,,^ F^(m, n) and compare it with the 1 /2 lower bound. Finally, 
it is instructive to extend our study to the case of using a different degenerate code to distill another 
type of entangled state subjected to another noise model, such as the Pauli channel [|16|l . 
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APPENDIX A: PROOF OF EQ. <HD 



We prove the validity of Eq. (|5D) by following the analysis in Sec. IIIIBI (And we follow the 
same notations as used in Sec. IIIIBl after possibly some straight-forward extension to the case of 
qudits.) First, we extend the definition of depolarization weight as follows. Let (fij, ary)"~g e K mn be 
a ordered rc-tuple. Then its depolarization weight is defined as the Hamming weight by regarding 
this n-tuple as a vector of elements in K. Clearly, Pr((5, 7) A s) = Pr((£(s, 5, 7)) where 



n-l 



6, 7) = {(J3j, ajfrl e K™ : 6 = £/3;, 7 = <*o, <*t = s { + a for t = 1,2, . . . ,n - 1} . (Al) 

Using the same argument as in the case of q = 2, we conclude that for k > 0, 

k\ 



1 

w{&{s,S,y);x,y) = - ^ 



q [a,,,) PI ai ' tl [\ ai '° l 

ieK, (gK* 

^-V[( ? -l)x + yf-* 



(n - kV 



n-k-a $ yi-a ,o^ao,o 



(A2) 



For A; = 0, we have to use a slightly different method to compute the depolarization weight enu- 
merator. By substituting Xq = y and x n = x for all rj ^ into the identity 



we have 



^ w^s*, 5,7); x,y) = [(q - \)x + y] n 



(A3) 



(A4) 



If K = GF(q), then by putting x = y and x v = (Op P x for all 77 ^ where p e GF(q)* into 
Eq. (IA3I) . we have 

2 o^V (£(£ 6, 7); x, y) = (y - xf . (A5) 

SeGF(q) 

If K = Z/qZ, then we put jco = y and x n = (Jlfx for all 77 ^ where p 6 (Z/qZ)* into Eq. (IA3I) . we 
arrive at 

£ <w((E(£ 5, 7); x, y) = (y - x)' 1 . (A6) 

SsZ/gZ 
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Consequently, we conclude that for K = GF(q) or Z/qZ, 

(q k -^[(q-l)x + y] n - k 



w(<£(f,d,yy,x,y) 



if k > 0, 



-{[(q-l)x + y] n + (q-l)(y-x) n } if k = and 8 = 0, (A7) 



-{[(^- l)jc + yf-Cy-jc)"} 



if k = and 5*0. 



Surely, for depolarizing channel, Pr((£, y) A it) = vv((£(s*, 5, 7); (1 - F)/(q m - 1), F) and 



Pr(s) = £ 







<f ' - 1 





F + 



(q-l)(l-F) 



(q m ~ 1) 



q m - 1 



1=0 



where 



/Ki) = e K™- 1 :k(lt) = i} 
Combining Eqs. (fATt -dA^t, we arrive at Eq. (I5TT) . 



(A8) 

(A9) 
□ 
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